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Abstract 


In gauge-mediated supersymmetry breaking theories the next-to-lightest su¬ 
persymmetric particle can decay during or after the nucleosynthesis epoch. 
The decay products such as photons and hadrons can destroy the light el¬ 
ement abundances. Restricting the damage that these decays can do leads 
to constraints on the abundance and lifetime of the NLSP. We compute the 
freezeout abundance of the NLSP by including all coannhilation thresholds 
which are particularly important in the case in which the NLSP is the lightest 
stau. We find that the upper bound on the messenger scale can be as stringent 
as 10 12 GeV when the NLSP is the lightest neutralino and 10 13 GeV when the 
NLSP is the lightest stau. Our findings disfavour models of gauge mediation 
where the messenger scale is close to the GUT scale or results from balancing 
renormalisable interactions with non-renormalisable operators at the Planck 
scale. When combined with the requirement of no gravitino overabundance, 
our bound implies that the reheating temperature after inflation must be less 
than 10 7 GeV. 
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1. Gauge-mediated supersymmetry breaking theories provide an interesting alternative 
to the usual gravity-mediated scenarios in transmitting supersymmetry breaking effects to 
the low energy world [IJ (for a review, see Ref. |]2|). Among the most attractive features of 
these theories is the natural suppression of flavour-violating interactions which is guaranteed 
by the gauge symmetry. In addition low-energy phenomenology is now governed by the fact 
that the gravitino is the lightest supersymmetric particle and this may lead to interesting 
collider signals. 

The next-to-lightest supersymmetric particle (NLSP) also plays an important role in 
low-energy phenomenology. Since the soft mass of a supersymmetric particle is determined 
by its gauge quantum numbers, the NLSP will be either a neutralino N { or the lightest, 
mainly right-handed stau fi, depending on the choice of parameters (the possibility of a 
sneutrino NLSP is now marginal). The relevant parameters which define the gauge-mediated 
model are the messenger index N (twice the sum of the Dynkin indices of the messenger 
gauge representations), and the supersymmetric and supersymmetry-breaking messenger 
mass parameters M and F. We also distinguish the supersymmetry breaking scale F felt 
by the messenger fields and the fundamental scale of supersymmetry breaking F 0 which 
determines the gravitino mass and couplings by defining k = F/F 0 < 1. In our analysis, we 
will trade F for the NLSP mass which is a parameter of more direct physical meaning. The 
Higgs mass parameters /i and B M are new independent inputs, which can be determined by 
imposing electroweak symmetry breaking, and therefore be fixed in terms of tan /3 and the 
algebraic sign of fi. As a result, the lightest neutralino turns out to be mainly R-ino. More 
details on the definition of the parameters can be found in Ref. |2]. 

The NLSP is not stable and eventually will decay into the gravitino. If these decays 
occur during the nucleosynthesis epoch the light element abundances can be drastically 
altered |3| For example, if the NLSP is a neutralino the dominant decay mode produces 
a photon which affects the primordial 4 He abundance. On the other hand hadronic decays 
can also prove dangerous for the nucleosynthesis products. In the case of the neutralino 
NLSP, the photon can hadronise or, in the case of a stau NLSP, the semileptonic decay 
of a tan lepton can produce hadronic showers which leads to energetic nucleons. Even 
though nucleosynthesis may be over, these energetic nucleons destroy 4 He or synthesise 3 He 
and tritium leading to the overproduction of the light elements. Alternatively if the NLSP 
actually decays hadronically during nucleosynthesis these nucleons will instead establish 
thermal equilibrium with the surrounding plasma by colliding with the ambient protons and 
4 He leading to the eventual increase of the n/p ratio. This results in a greater abundance 
of 4 He, 3 He and deuterium D. 

In order to avoid the destructive effects on the nucleosynthesis products the lifetime of 
the NLSP must be restricted so as to decay sufficiently well before it can interfere with the 
nucleosynthesis products or if it decays during nucleosynthesis that the enhanced light ele¬ 
ment abundances are consistent with astrophysical observations. In addition the abundance 
of the NLSP at the time of decay will also be important. Since the abundance and lifetime 
are related to the messenger scale M, an upper bound can be placed on the messenger scale 
(for a typical set of the other gauge-mediated parameters), depending on whether the neu¬ 
tralino or stau is the NLSP. These bounds will be shown to be fairly generic for most of the 
parameter space of gauge-mediated theories. 

2. The damaging effects of the NLSP decay products during the nucleosynthesis epoch 
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constrains the abundance and lifetime of the decaying particle. In order to obtain a bound 
on the NLSP lifetime a detailed calculation of the NLSP abundance at the time of decay 
must be performed. This amounts to calculating the NLSP abundance at the time of freeze 
out when the NLSP is no longer in chemical equilibrium. 

We will consider the two separate cases of a neutralino and stau NLSP with mass utnlsp- 
For moderate values of tan /?, f\ is lighter than the neutralino whenever 
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For large tan (3 this region becomes slightly larger because of the stau left-right mixing. In 
order to determine the NLSP abundance we consider all relevant channels which change the 
number of NLSP’s. When the NLSP is a neutralino the most relevant annihilation channels 
are those consisting of fermions and gauge bosons as depicted in Table |. The annihilation 
channels into Higgs bosons are suppressed because the lightest neutralino is mainly H-ino. 
For the same reason, the neutralino NLSP generically is not degenerate in mass with other 
particles, so it will not be important to consider coannihilations for this case. We can also 
neglect annihilation channels which involve gravitino vertices since these lead to scattering 
amplitudes which are suppressed by a factor i r n‘^ LSP /F 0 . Since we are interested in NLSP 
lifetimes of the order of the nucleosynthesis timescale, this suppression makes all gravitino- 
emission processes negligible. 

On the other hand when the NLSP is the stau, the lightest smuon and selectron have a 
small mass difference with the NLSP: 


ra A i,ei - rn fl ^ m 2 

Wr, 2m~ 

Here A r is the supersymmetry-breaking trilinear term and r are the flavour-independent 
contributions to the left and right slepton masses (including the .D-term contribution). Since 
in gauge-mediated theories p 2 > m~ , the mass difference in Eq. (§) is always positive. 
Because of the approximate mass degeneracy among the sleptons, the calculation of the 
stau relic abundance must include all coannihilation processes listed in Table |j. We also 
need to compute the smuon and selectron density at the decoupling time, since these “co- 
NLSP’s” are also responsible for producing damaging decay products. 

The NLSP abundance is determined by considering the evolution of the number density n t 
of particle i which is governed by the Boltzmann equation. In the presence of coannihilations 
the Boltzmann equation can be written as [jllj 

= -3 Hn - ( a eS v ) (n 2 - n 2 q ) (3) 

where H is the Hubble expansion parameter, n = YU n i an d the thermal average of the 
effective cross section is defined as 
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TABLE I. Final state annihilation channels including coannihilations in parentheses. 


initial state 

final state 

JV 1 JV 1 

zz,w+w-Jf 


ZZ, W+W- , 77, Z 7 , hh^h, Zh , ff 

U ±f i 

t ± t ± , e ± e ± ) 


(r ± g ± ,T ± e ± ,g ± e ± ) 

(fifii > if&i . hfef ) 

(r ± /r =F , VrVp), (r ± e =F , n T n e , n T n e ), (p ± e =F , u^v e ) 


The individual cross sections a t] include the processes listed in Table [I]. v l3 is the relative 
velocity and n eq is the equilibrium number density. 

The cross sections for all annihilation channels are numerically evaluated using the 
CompHEP software package [[□]]. After performing the thermal average and including all 


relevant coannihilation thresholds, the ratio U eq of the equilibrium number density n eq to 
the entropy density s is given by 
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where x = m^sp/T, gi is the number of internal degrees of freedom and g e g = 81. The 
NLSP abundance Tnlsp at the time of nucleosynthesis is given by Lnlsp = Y eq (T F ), where T F 
is the freezeout temperature. The freezeout temperature is determined from the condition 
n eq (v e ftv) = H(T). 

Thus in the case of the neutralino NLSP the abundance Inlsp is just simply while 
for the stau NLSP it will be Inlsp = If + Yp : + Yg which is the sum of all the “co-NLSP” 
abundances. The effect of including the coannihilating channels between the stau and the 
“co-NLSP’s” changes the pure stau abundance by up to ~ 50%. This can be seen in Fig. [I] 
where the stau abundance with and without coannihilations is shown. Notice the important 
reduction of the relic abundance for values of close to half the Higgs mass, caused by 
the resonant annihilation channel. It should be noted, however that given the present LEP 
bound on stable fi, this possibility is no longer realistic for the stau. 

3. Let us now consider the effect of the NLSP decaying during the nucleosynthesis epoch. 
Suppose first that the NLSP is a neutralino. The dominant decay mode of the neutralino is 
into a photon and gravitino with a decay rate given by H JT^Jl3|] 
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where cos Qw + Ai 2 sin 9w \ 2 and AY 12 are the NLSP gaugino components in stan¬ 

dard notation. There are also other decay modes which involve the Z or a Higgs boson, but 
these modes are suppressed by a f3 8 phase-space factor. The photon of the dominant decay 
mode has a dramatic effect on the light element abundance 1 . I 11 the radiation dominated 


1 Notice that the gravitinos released in the NLSP decays do not thermalise and their present 
contribution to the energy density of the Universe is negligible. 
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FIG. 1. Comparison of the stau abundance with and without coannihilation (lower and upper 
curves, respectively). The gauge-mediated parameters are M = 10 13 GeV, N = 12, tan/3 =1.1 
and sgnp = 1 . 


thermal background high energy photons initiate electromagnetic cascades and create many 
low-energy photons capable of photodissociating the light elements 0. These photodisso¬ 
ciation effects caused by the photon of the decaying neutralino can be parametrised by a 
“damage” factor 


dJ = TTTTVi 
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where 77 is the baryon-to-photon number density ratio. This quantity is constrained 
by astrophysical observations of light elements. Recently new observations of the D and 
4 He primordial abundance have led to a reanalysis of the constraints on the abundance 
and lifetime of long-lived particles jl|. In particular, controversy in the measurement of 
primordial deuterium have led to two different values for the deuterium abundance Xp 
(normalised to hydrogen). In Ref. [17I| a high value of the deuterium abundance Xp = 
(1.9 ± 0.5) x 10 -4 is quoted while in Ref. |R| a low value X D = (3.39 ± 0.25) x 10 


-5 


is 


obtained. 

The neutralino abundance Y Nl depends on the details of the sparticle spectrum and can 
be calculated for a particular choice of the parameters N, M , F, tan /3 and sgn/ 1 . Using the 
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constraints in Ref. || we show in Fig. |2jthe bound on the messenger scale for a representative 
neutralino NLSP scenario in which N — 2, tan/3 = 2 and sgn/i = 1, for both high and low 
values of the measured X D . We have also assumed that k = 1 and the limiting value of M 
grows approximately linearly with k. Notice that the high deuterium value gives a slightly 
more stringent bound than the low deuterium value, because for lifetimes less than 10 6 
seconds deuterium is effectively photodissociated. 



FIG. 2. Upper bound on the messenger scale M as a function of the neutralino NLSP mass 
mjVi from photodissociations where k = 1, N = 2, tan/3 = 2 and sgnp = 1. Both bounds use a 4 He 
abundance Yp( 4 He) = 0.234 but the dashed line uses a low value of the deuterium abundance while 
the solid line uses the high value of the deuterium abundance (see text). 

The bound will change slightly for different values of the gauge-mediated supersymmetry 
breaking parameters, but typically for the neutralino NLSP scenario lifetimes greater than 
10 5 seconds are ruled out by photodissociations. 

Note that photodissociations are only important after nucleosynthesis is over (which 
corresponds to times later than 10 4 seconds). This is because during earlier times when 
thermal photons are more energetic and numerous, photon-photon interactions are much 
more probable than photon-nucleus interactions. Consequently stronger constraints can 
only be obtained if we consider processes that affect nucleosynthesis at times earlier than 
10 4 seconds. 
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However before discussing the relevant processes for times earlier than 10 4 seconds it is 
possible to obtain stronger constraints than those from photodissociations by considering 
hadronic decays of the NLSP ||. If the NLSP decays hadronically during times > 10 4 
seconds then the photodissociation bound needs to be reconsidered since while D is destroyed 
by photodissociation it can be produced by the hadronic showers. In fact even if the NLSP 
decays exclusively into photons it is possible that hadronic showers will be generated. The 
main effect of hadronic showers is not only to increase the abundance of D but also the other 
light nuclei 3 He, 6 Li and 'Li. In particular D and 3 He arise mainly from “hadrodestruction” 
processes, which occur when an energetic nucleon breaks up the ambient 4 He nucleus into D 
and 3 He. The other light nuclei, b Li and 'Li arise from “hadrosynthesis” of 3 He, T or 4 He. 
More complete details can be found in Ref. ||. Again the “damage” of the hadronic decays 
on the primordial nuclear abundances can be parametrised as 


dg = d^B 


7 ■ LJ h 
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where = (i'B/5)BhJ r is an effective baryonic branching ratio which depends on the true 
baryonic branching ratio Bh , the baryonic multiplicity z/g and a factor T representing the 
dependence of the yields on the energy of the primary shower baryons. 

One of the contributions to the neutralino hadronic branching ratio comes from the decay 
Ni —> ZG 


B h ~ 0.7 tan 2 9 W 1 
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where the hadronic branching ratio of the Z-boson is 0.7, and we have assumed a pure R-ino. 
Although this contribution, and the analogous one from Higgs decay, can be forbidden by 
phase space, there is always at least a contribution to Bh of order a from photon conversion 
into a qq pair. The explicit expression for Bh can be obtained from Eq. (13) of Ref. [IB . 
Using the corresponding value of dg, one Ends that the overproduction of 'Li constrains the 
lifetime of the neutralino to be shorter than 10 4 seconds. This is also true if the NLSP is 
the stau. The stau decays predominantly into a tau and gravitino and since the tau has a 
large hadronic branching ratio, the corresponding value of dg is much larger than that for 
the neutralino. Consequently the lifetime of the stau must also be less than 10 4 seconds if 
the overproduction of 'Li is to be avoided ||. 

It is clear that in order to obtain precise constraints on the abundance and lifetime of 
the NLSP we need to consider processes occuring at times earlier than 10 4 seconds. At these 
times the main decay products that interfere with nucleosynthesis are hadrons. Hadronic 
showers induce interconversions between the ambient protons and neutrons thus changing 
the equilibrium n/p ratio [Q. In particular during the lifetime interval Tnlsp ~ 1 — 100 
seconds the overall effect of the hadronic decays is to convert protons into neutrons. The 
additional neutrons that are produced are all synthesised into 4 He and thus hadronic decays 
during this time interval are constrained by the observational upper bound on the primordial 
helium abundance Tp( 4 He). 

Eventually the neutron fraction falls to zero because all neutrons are contained in the 
4 He nuclei and the remaining neutrons created by the NLSP decay increase the deuterium D 
abundance. Furthermore for tnlsp ~ 100 — 1000 seconds the 3 He abundance is also increased 
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by D-D burning. After Tnlsp ~ 10 4 seconds all the neutrons arising from NLSP decay will 
themselves decay before forming D. Thus in the interval 10 2 — 10 4 seconds the appropriate 
constraint on hadronic decays arises from the observational bounds on (D + 3 He)/H. 

The overall effect of these hadronic decays has been considered in Ref. [f||, where the 
constraints on the abundance and lifetime of the decaying particle are parametrised by 

r _ NjetBfr ( n{Ej e t )) , . 

J ~ 2 (n(33 GeV)) 1 J 

where Nj et is the number of jets, Bh is the hadronic branching ratio and ( n(Ej et )) is the 
average charge multiplicity for a jet with energy Ej et and is given by 


(■ n(Ej et )) = 1 + 0.0135 exp 
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Since at the parton level the neutralino NLSP decays into three particles we will assume 
that Ej et = mNLSp/3 and N jet is the number of quarks at the parton level. In the case of 
the stau in which the tau decays semi-leptonically we assume Ej et = mNLSp/6- 

The most stringent constraints on the NLSP abundance and lifetime in the interval 
10 2 — 10 4 seconds come from the primordial deuterium abundance. Previous constraints 
on Lnlsp/ an d the lifetime were obtained for X D < 10~ 4 and rj — 3 x 10~ 10 ,1CT 9 0. In 
order to use the more recent measurements, we rescale the previous constraints for the new 
values X D = (1.9 ± 0.5) x 10 -4 or X D = (3.39 ± 0.25) x 1CT 5 . This rescaling can be done 
analytically by using the change in the deuterium abundance J3[] 


XX D = 


/^Inlsp 
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where Y H is the hydrogen density, ep is the fraction of injected neutrons that end up in 
deuterium and a n is the number of nn pairs per NLSP decay. This equation is valid as long 
as Xo -C e d, which conveniently holds whenever the deuterium constraint is important. 
Since the limit comes from the overproduction of deuterium and the standard Xd prediction 
decreases with , we make the most conservative choice of rj — 6 x 10 _1 °, which is the largest 
value compatible with 4 He and °Li abundances. 

When the high value of the deuterium abundance is used, one obtains no constraints in 
the region 10 2 — 10 4 seconds. All values of the NLSP abundance and lifetime are consistent 
with high X D value and consequently there is no improvement on the lifetime upper bound 
of 10 4 seconds obtained from the overabundance of 'Li. 

On the other hand stringent constraints in the interval 10 2 — 10 4 seconds are obtained 
when the low value of the deuterium abundance is used. Let us consider the two NLSP 
scenarios separately. First, when the NLSP is the neutralino the scaling parameter / is 
determined using Nj et = 2 and B h = 1CT 2 . Rescaling the solid curve of Fig. 4 in Ref. Q 
enables one to determine the constraint arising from the low value of X D and ij = 6 x 10~ 10 . 
Thus calculating the value for the neutralino number density for a generic set of gauge- 
mediated supersymmetry breaking parameters leads to constraints on the abundance and 
lifetime of the neutralino which can be expressed as a bound on the messenger scale M, see 
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FIG. 3. Upper bound on the messenger scale M as a function of the neutralino NLSP mass 
from hadronic decays where k = 1, N = 2, tan/3 = 2 and sgnp, = 1. The solid (dashed) 
line corresponds to the bound assuming a low (high) deuterium measurement. More specifically the 
dashed line represents the 10 4 second lifetime contour which arises from the 7 Li overabundance. 


Fig. |3j. The bound on M does not significantly depend on the value of tan (3, but grows 
approximately linearly with N. 

Let us now discuss the case of the stau NLSP. At the decoupling time, the abundances 
of smuons and selectrons are comparable to the stau abundance. The cosmological fate of 
the frozen-out smuons and selectrons depends on the mass difference in Eq. (|2|). If tan (3 is 
large, then mi — m fl > m T — (i = e,/i) and the decay — > f\rf is kinematically open. 
The value of tan (3 for which this transition occurs depends on the parameter choice, but it 
is typically of order 10, as it can be estimated from Eq. (]2|). If this mode is accessible, it 
dominates over the two-body decay into gravitino. This can be simply understood from the 
decay rate expression [|T7[] in the limit of vanishing lepton masses and small slepton mass 
difference, 
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Here m^ is the mass of the H-ino which mediates the decay, and is assumed to be much 
larger than mg . In this case, smuons and selectrons decay soon after decoupling, with 
each one producing a stau in the final state. Therefore, the relevant stau abundance at the 
nucleosynthesis epoch is determined by the sum over the three slepton abundances at the 
decoupling time. The stau hadronic branching fraction comes from the semileptonic tau 
decay, and it is given by B h ~ 0.65. 

For smaller values of tan f3 the neutralino-mediated three-body process is forbidden, and 
the rates of the two competing slepton decay modes are 

k 2 rrfe 

r - eG) = (15) 


C 2 M ' 4 . 

T(^ -»■ rf vgpr) ~ - m fl ) 5 sin 2 9 f sin 2 9g. (16) 

Iotz 6 mz+ 1 
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Here m x + is the chargino (gaugino) mass which mediates the decay, and 9f, 9g are the 
left-right mixing angles in the slepton system. The process in Eq. ([15]) is suppressed by 
the large value of F required in our study of nucleosynthesis and the process in Eq. (|T6| ) 
is suppressed by the small mixing angles proportional to the corresponding lepton masses. 
It turns out that in the small tan/? regime and for mg ~ 100 GeV, the two rates are 
comparable when \[F ~ 10 9 GeV, which is just the value of \f~F necessary to have a stau 
decay during nucleosynthesis. Therefore, depending on the particular choice of the gauge 
mediation parameters, we can encounter different situations. The first possibility is that 
T(G —> IG) > T{t[ —> Tfz^P r ), in which all light sleptons decay directly into gravitinos 
around the same time. In this case smuons and selectrons do not significantly contribute 
to the effective hadronic branching fraction, because real electrons and muons cannot decay 
into hadrons. As tan f3 is increased and mg is decreased, we go first to a regime in which 
r(/H H _z/ /i^r) > r(G —> £G) > T(eq —> ff/y e P T ) and then in a regime in which the 
three-body decays dominate for both fi\ and e\. The two regimes are possible because the 
decay process in Eq. (|l(|) depends on the lepton mass through the left-right mixing angle. 
In the first regime, only the smuon contributes to the effective hadronic branching fraction, 
while in the second one both the smuon and selectron contribute. 

With respect to a neutralino NLSP of equal mass, a stau NLSP has a larger hadronic 
branching ratio Bh , but a smaller relic abundance because of the additional annihilation 
channels. The two effects roughly compensate each other, although the bound on the NLSP 
lifetime is slighly weaker in the fi case. When expressed in terms of M, the bound appears 
even weaker because we have to choose a very large value of N in order to satisfy the 
requirement of Eq. (|f[) for a stau NLSP. The bound on the messenger mass for N = 12 and 
tan/3 = 1.1 is shown in Fig. f|. For this choice of parameters, the selectron and the smuon 
dominantly decay directly into gravitinos. In the figure we only show the bound arising from 
the low deuterium value. The bound for the high deuterium value corresponds to a lifetime 
which cannot be achieved for this choice of N. 
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FIG. 4. Upper bound on the messenger scale M as a function of the stau NLSP mass mf from 
hadronic decays where k = 1, N = 12, tan/3 =1.1 and sgn/j = 1. 


4. Although the mass spectrum of supersymmetric theories with gauge mediation can 
be predicted in terms of a few parameters, the experimental and cosmological features of the 
theory can drastically change as the messenger mass M is varied in the allowed range between 
about 100 TeV and 10 15 GeV. It is thus very important to study possible constraints on the 
value of M. The requirement that the NLSP decay products do not upset the successful 
predictions of standard big-bang nucleosynthesis provides such a constraint. Here we have 
studied in detail the effects of the decaying NLSP on the light-element primordial abundances 
and shown that the most dangerous damages come from the hadronic decay modes. After 
computing the NLSP relic abundance, including coannihilation effects, we conclude that the 
injection of hadronic jets in the Universe at times later than 10 4 seconds grossly overproduces 
'Li. This is true for both neutralino and stau NLSP. At earlier times, the most relevant 
bounds come from deuterium overproduction. At present, there is a disagreement on the 
extracted observational value of the primordial deuterium. Using Xp = (3.39 ± 0.25) x 
10 -5 |n|, no further limit on the NLSP is derived, while with Xo = (1.9 ± 0.5) x 10 -4 [L4[] 
the limit from 'Li can be improved. The corresponding limits on the messenger mass M 
are shown in Figs. [3] and [| for a representative choice of gauge mediation parameters, in the 
case of neutralino and stau NLSP respectively. 
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It should be noticed that the bounds on the messenger mass presented here can be evaded 
in the presence of other interactions which can lead to a fast NLSP decay, e.g. -R-parity 
violation. Nevertheless these bounds have significant implications for model building. In 
particular they disfavour models in which M is close to the GUT scale or models in which the 
messenger scale results from balancing renormalisable interactions with non-renormalisable 
operators at the Planck scale HT8|,|T9i . 

It is also interesting to notice that the nucleosynthesis bound discussed here is comple¬ 
mentary to the bound obtained from gravitino overabundance. Indeed when mg, is larger 
than about a keV, gravitinos which were in thermal equilibrium at early times, give a con¬ 
tribution to the present energy density larger than the critical value. It is then necessary to 
assume that gravitinos have been diluted by some mechanism. Let T max be the temperature 
at which the ordinary radiation-dominated Universe begins. This corresponds to either the 
reheating temperature after an inflationary epoch or to the temperature associated with sig¬ 
nificant entropy production. The requirement that gravitinos do not overdose the Universe 
gives the following contraints on T max [|20|,pT|| 

Tmax ~ 100 GeV — 1 TeV for 2 h 2 keV ^ mg ^ 100 keV, (17) 


/ \ /TeV\ 

T max ^ 10 TeV h 2 ( G ) - for mg ^ 100 keV, (18) 

V100 keV / \rrig J 

where h is the Hubble constant in units of 100 km sec -1 Mpc -1 and rrig is the gluino mass. 
To compare this bound with the nucleosynthesis result, it is convenient to express it in terms 
of the messenger mass M, the messenger index N, and the pure H-ino mass m ^. We find 


- max rsj 


< 100 GeV - 1 TeV for 10 8 h 2 GeV < 
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m N , 


kN V100 GeV 


) < 5 x 10 9 GeV, (19) 


T max £ 5 X lO" 6 h 2 — 
kN 



for — 
kN 


m Nl 


100 GeV 


> 5 x 10 9 GeV. 


( 20 ) 


This shows that there is no gravitino problem as long as 
M ^ 10 8 h 2 GeV fclV(100 GeV/mjVj)- For larger values of M, there is a very stringent 
constraint on T max which requires inflation at particularly low temperatures. As M grows 
this limit becomes weaker, but eventually the nucleosynthesis bound on the messenger mass 
sets in. When the two bounds are combined, we find that the reheating temperature after 
inflation should be less than about 10' h 2 GeV which is stronger than the bound ~ 10 10 
GeV usually obtained in gravity-mediated scenarios |5|j|j2(j. Our limit is valid for gauge- 


mediated theories in which the LSP gravitino is heavier than a keV and in which there are no 
new interactions uncorrelated with the supersymmetry-breaking scale mediating the NLSP 
decay (like A-parity violating interactions). 
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